Astron. Nachr. / AN ?, No. ?, 0-[II](??) / DOI please setDOI! 



o 

(N 



o 

u 



>< 



X-ray luminosity functions of different morphological and X-ray type 
AGN populations 

M. Povic^^,^ A. M. Perez Garci'a^^, M. Sanchez-Portal^^, A. Bongiovanni^^^, J. Cepa'' \ M. 
Fernandez Lorenzo^^, M. A. Lara-Lopez^ ^, J. Gallego*^, A. Ederoclite^^, I. Marquez^, J. Masegosa^, 
E. Alfaro^, H. Castaiieda^, J. I. Gonzalez-Serrano^", and J. J. Gonzalez^^ 



CO , ^ Instituto de Astroflsica de Canarias (lAC), La Laguna, Tenerife, Spain 

Institute de Astroflsica de Andalucfa (lAA-CSIC), Granada, Spain 
Asociacion ASPID, Apartado de Correos 412, La Laguna, Tenerife, Spain 
Departamento de Astroflsica, Universidad de La Laguna (ULL), La Laguna, Spain 
' ^ Herschel Science Centre, ESAC/INSA, Madrid, Spain 

D : Australian Astronomical Observatory (AAO), Sidney, Australia 

Departamento de Astroflsica y CC. de la Atmosfera, Universidad Complutense de Madrid, Madrid, Spain 
' * Centro de Estudios de Fsica del Cosmos de Aragon (CEFCA), Teruel, Spain 

^SJ ■ ^ Escuela Superior de Fisica y Matematicas, Intituto Politecnico Nacional (IPN), Mexico D.F, Mexico 

Instituto de Fisica de Cantabria, CSIC-Universidad de Cantabria, Santander, Spain 
Instituto de Astronomfa UNAM, Mexico D.F, Mexico 



10 

11 



Received ?, accepted 



1 



(~| ' Published online later 
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^2 ' Luminosity functions are one of the most important observational clues when studying galaxy evolution over cosmic time. 

t/^ ' In this paper we present the X-ray luminosity functions for X-ray detected AGN in the SXDS and GWS fields. The 

limiting fluxes of our samples are 9.0 x 10^^^ and 4.8 x 10~^® ergcm"^ s^^ in the 0.5-7.0keV band in the two fields, 
respectively. We carried out analysis in three X-ray bands and in two redshift intervals up to z < 1.4. Moreover, we derive 
the luminosity functions for different optical morphologies and X-ray types. We confirm strong luminosity evolution in all 
^ , three bands, finding the most luminous objects at higher redshift. However, no signs of density evolution are found in any 

. tested X-ray band. We obtain similar results for compact and early-type objects. Finally, we observe the 'Steffen effect', 

where X-ray type-1 sources are more numerous at higher luminosities in comparison with type-2 sources. 

^ © ?? WILEY-VCH Verlag GmbH & Co. KGaA, Wcinhcim 
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' 1 Introduction the availability of data for large samples of sources, cov- 

ering a wide range of redshift and luminosity, in order to 
The luminosity function (LF) is among the fundamental accurately determine the shape of the LF and its evolution 
tools for studying the formation and evolution of galaxies cosmic time. 



(e.g. Schechter 1976). Galaxies with an Active Galactic Nu- X-ray emission is one of the principal characteristics of 

y ; cleus (AGN) in their centres are among the most luminous AGN activity. Since X-ray data are relatively unaffected by 

- -' sources in the Universe. Therefore, they provide unique in- absorption, they have been shown to be very useful for se- 

formation about the properties of high-redshift galaxies that lecting large AGN samples. Hard and very-hard X-ray sur- 

is important for understanding how galaxies formed and veys are especially important, since they are able to detect 

evolved. However, the manner in which a galaxy enters the significant populations of sources that are obscured in opti- 

active phase, the duration of the activity, and the origin of cal and NIR bands. Therefore, X-ray data provide the most 

the accreting material have not been clarified yet. The use complete and unbiased samples of AGN, and the X-ray lu- 

of LPs to determine the distribution and evolution of the minosity function (XLF) can provide important knowledge 

AGN accretion history has been shown to be essential for about AGN evolution over cosmic time (e.g., Barger et al. 

constraining models of formation and evolution of super- 2005; Barger & Hasingeret al. 2005). 

massive black holes (SMBHs), AGN feedback, fuelhng me- The AGN LF has been widely studied over past years, 

chanisms, and the relation between active and non-active providing important constraints on models of galaxy for- 

galaxies (e.g., Aird et al. 2010, La Franca et al. 2010, Fiore mation and evolution. In the last decade, deep X-ray (e.g., 

et al. 2012 and references therein). LF computation requires Ueda et al. 2003; La Franca et al. 2005; Hopkins, Richards, 

& Hemquist 2007; Silverman et al. 2008a, 2008b; Ebrero et 

* Corresponding author: e-mail: mpovic@iaa.es al. 2009; Yencho et al. 2009; Aird et al. 2010 and references 

therein), optical (e.g., Boyle et al. 2000; Fan et al. 2001, 



®WILEY 

) interScience' 



© ?? WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim 



Astron. Nachr. / AN (??) 



1 



2004; Croom et al. 2004; Richai'ds et al. 2005, 2006; Jiang 
et al. 2006; Cool et al. 2012), IR (Dai et al. 2009 and ref- 
erences therein; Wu et al. 201 1 and references therein) and 
radio (e.g.. Kaiser & Best 2007; Kellermann et al. 2012 and 
references therein) surveys studied large samples of AGN 
over a broad range of redshift (z < 6). These multiwave- 
length studies have provided a wider view of AGN prop- 
erties. As akeady mentioned. X-ray selection provides one 
of the most complete samples of active galaxies; however 
they may miss the most obscured (Compton thick) sources. 
These sources can be observed in the IR. In contrast, optical 
data are important for providing the redshift and AGN mor- 
phology. Current observations show that the 
AGN space density reaches a broad peak at z ?» 1 - 3 (e.g. 
Hasinger et al. 2005; Silverman et al. 2008b; Aird et al. 
2008; Brusa et al. 2009) followed by a decUne of w2dex 
to the present day (e.g., Ueda et al. 2003; Hasinger et al. 
2005; Barger et al. 2005). Moreover, X-ray studies (e.g.. 
Page et al. 1997; Miyaji et al. 2000, 2001; La Franca et 
al. 2002; Cowie et al. 2003; Ueda et al. 2003; Hasinger 
et al. 2005; Aird et al. 2010 and references therein), but 
also optical, IR, and radio measurements (e.g.. Hunt et al. 
2004; Cirasuolo, Magliocchetti, & Celotti 2005; Matute et 
al. 2006) indicate that the space density of low-luminosity 
AGN peaks at lower redshift than the bright AGN. In order 
to explain this effect, called 'cosmic downsizing' (Cowie et 
al. 2003; Merloni 2004; Heckman et al. 2004, Babic et al. 
2007; Fanidakis et al. 2012), most of the proposed models 
of AGN formation and evolution consider different forms of 
AGN feedback processes as the responsible mechanism (Di 
Matteo et al. 2005; Hopkins et al. 2006; Croton et al. 2006; 
Hopkins et al. 2008). 

Many studies of the AGN LF have been published, lead- 
ing to different models of AGN evolution. However, their 
analyses were restricted to the evolution of the entire AGN 
population. In this work, in addition to a study of the en- 
tire AGN sample, the X-ray luminosity function was also 
studied for different morphologies and X-ray obscurations 
for redshifts z < 1 .4. The results presented here represent a 
pilot study of a much larger sample of AGN that will be car- 
ried out in the COSMOS field (Scoville et al. 2007) with the 
final aim of studying the evolution of different AGN trigger- 
ing mechanisms over cosmic time. The paper is structured 
as follows: in Section |2l we describe the data used in this 
work, including X-ray (Section [2. It , optical (Section [2. 2| i, 
near- and mid-infrared (NIR and MIR) data (Section |2.3l l, 
and provide brief summary of photometric redshift and k- 
correction estimations of the optical counterparts of X-ray 
emitters (Section lZ^l i. In Section[3] we describe the optical 
morphology and X-ray type classification of X-ray selected 
AGN in both fields. In SectionlH we present the analysis of 
the XLFs, with the luminosity function measurements de- 
scribed in Section|4T|and discussion of the LFs of the whole 
AGN sample in Section |42] of the different morphological 
types in Section |431 and of the different X-ray types in Sec- 
tion |4.4| Our study of the evolution of X-ray type- 1 sour-ces 



is given in Section |431 Finally, the results obtained in this 
work are summarized in Section |5] 

The concordance cosmology with S1a = 0.7, i7j\/=0.3, 
and Ho = 70 km s~^ Mpc~^ is assumed. Unless otherwise 
specified, all magnitudes are given in the AB system (Oke 
& Gunn 1983). 

2 Observational data 

This work involves sources belonging to two fields: Su- 
baru/XMM-A^evvfo« Deep Survey (SXDS) and Groth-West- 
phal Strip (GWS). These fields were selected for the OSI- 
RIS Tunable Emission Line Object (OTELO) survey (Cepa 
et al. 2008). This is a deep on-going survey of emission 
line sources that are being observed with tunable filters in 
the Optical System for Imaging and low Resolution Inte- 
grated Spectroscopy (OSIRIS) on the Gran Telescopio de 
Can anas (GTCQ) telescope. Results presented in this paper 
complement studies of the AGN/QSO LFs carried out in the 
OTELO survey also using OSIRIS data. 

The SXDS is a large survey covering > 1 square degree 
centred at RA = 02 ''18'" and DEC= -05°00' (Sekiguchi et 
al. in prep), and observed from X-ray to radio. 

The GWS field is a part of the Extended Groth Strip 
(EGS) centred at RA = 14'' 17™ and DEC = + 52°30' . Broad 
band optical BVRI observations were performed by the 
OTELO group (Cepa et al. 2008). The EGS was also ob- 
served by the All-Wavelength Extended Groth Strip Inter- 
national Survey (AEGISS) team, with the aim of obtaining 
deep images and spectroscopic data from X-ray to radio, 
covering an area of 0.5 - 1 square degree (Davis et al. 2007). 

We mainly use X-ray and optical data in this work. NIR 
and MIR data were used to complement the optical infor- 
mation for the estimation of photometric redshifts estima- 
tions (see Section 1241 1. The observational data are briefly 
described below. 

2.1 X-ray data 

The SXDS field was observed with XMM-Newton at 
seven pointings in the 0.2-lOkeV energy range (PI Michael 
G. Watson) and covering ~ 1.3 deg^. The central observa- 
tion is deepest with a nominal exposure time of ^ lOOksec. 
It is surrounded by six shallower observations each with an 
exposure of ^ 50ksec. 

The required European Photon Imaging Camera (EP- 
IcH) data were collected from the XMM-Newton v.5.0 sci- 
entific archive (XSA0) and were processed by our team us- 
ing standard procedures of the Science Analysis System 
(SAsEt) v7.1.2 package, and the latest relevant Current Cal- 
ibration Files (CCF). Source detection was based on the 

' http://www.gtc.iac.es/ 
^ http://aegis.ucolick.org/ 

^ http://xmm.esac.esa.int/extemaI/xmm_user_support/documentation 
/technical/EPIC/ 

http://xmm.esac.esa.int/xsa/ 
^ http://xmm.esac.esa.int/external/xmm_data_analaysis 
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edetect chain SAS procedure. In order to minimize spu- 
rious detections the hkeUhood threshold parametei|^ was set 
at L= 14. Thus the probability that the source is real ex- 
ceeds P = 0.99 which implies less than one fake source per 
instrument, per pointing, and per band. A detailed descrip- 
tion of X-ray data reduction and source detection is avail- 
able in Povic et al. (2012). We detected 1121 sources with 
fluxes greater than our flux detection limit of 9.0 x 10~^^ 
erg cm^^ s^^ in the 0.5 - 7.0keV band. 

The GWS field was observed with Chandra covering 
three continuous fields centred at the original HST GWS. 
All observations were carried out by the AEGIS collabora- 
tion using the Advanced CCD Imaging Spectrometer (AC- 
IS-fl) instrument, with a total exposure time of ~ 200 ksec 
in each field, covering an area of 0.08 deg^. 

All datasets were gathered from the Chandra Data Arc- 
hivH and were processed by our team using the Chandra 
Interactive Analysis of Observations (CIAOfl v4.0 pack- 
age and Calibration Data Base (CALDB) v3.4.0. Standard 
reduction procedures were applied to produce new level 2 
event files, restricted to the 0.5-8 keV range in order to avo- 
id high background spectral regions. Good time intervals 
were defined by means of celldetect and a threshold 
of S/N = 3. Absolute astrometry was improved, where the 
average positional accuracy improvement is 6.3% with re- 
spect to the original astrometry provided by the spacecraft 
attitude files, and the event files were merged. We applied 
the CIAO wavdetect Mexican-Hat wavelet source detec- 
tion program to all bands. We set a significance threshold 
of 2 X 10^^, expecting about 12 spurious sources in the en- 
tire catalogue. The final catalogue has 639 unique X-ray 
emitters, with fluxes greater than the flux detection limit of 
4.8 X 10^1^ erg cm^^ s^^ in the 0.5-7.0keV band. A de- 
tailed description of X-ray data reduction and source detec- 
tion can be found in Povic et al. (2009). 

In this paper, we derive the XLF in the 0.5-2 keV (soft) 
and 0.5-7 keV (total) bands for both the SXDS and GWS 
samples, and in the 2-lOkeV (hard) band for the SXDS 
sample. To separate between the X-ray unobscured and ob- 
scured sources (see Section |3]l, we use HR(2 — 4.5keV/ 
0.5 — 2keV) hardness ratio defined as: 

^ CR{2-A.bkeV)-CR{Q.b-2keV) 
^ CR{2 - 4.5 keV) + CR{0.5 - 2 keV) ' 

where CR is the count rate in a given energy band. 

We tested the reliability of our X-ray number counts 
in the SXDS and GWS fields by comparing with number 
counts obtained in deeper (e.g. Chandra Deep Field North 

* L = -ln(l - P), where P is the probabihty to have a spurious detection 
due to random Poisson fluctuations. 

' http://cxc.harvard.edu/caI/Acis/ 

* http://asc.harvard.edu/cda/ 
' http://cxc.harvard.edu/ciao/ 
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Fig. 1 Cumulative logN-logS functions for the SXDS and 
GWS fields in the soft band, for all detected X-ray sources 
(red and blue thick solid lines, respectively) and for X-ray 
emitters with optical counterpart (red and blue thin dashed 
lines, respectively). Three distributions, with different ex- 
posure times and effective areas, are represented for com- 
parison: GWS by Nandra et al. (2005) (black dotted line), 
CDFN (black dash-dot-dashed line), and ELAIS-Nl (black 
dashed line). 

(CDFN); Alexander et al. 2003) and less deep (e.g., ELAIS- 
Nl; Manners et al. 2003) surveys. Figure [T| shows the cu- 
mulative number count distributions (logN-logS) per square 
degree in the soft band. Two curves have been computed 
for the SXDS and GWS sample, one for the complete sam- 
ple of X-ray detections and the other for the sample of X- 
ray detections with optical counterparts. For objects brighter 
than ~4 x lO^^'^ erg s^^ cm^^ in the case of SXDS, and 
^ 10~^^ergs~^ cm~^ in the case of GWS, both distribu- 
tions (complete X-ray sample and X-ray sample with opti- 
cal counterparts) coincide, while for fainter objects the den- 
sity of X-ray sources with optical counterparts starts to de- 
crease. In contrast, the difference in the density of X-ray 
sources between the distributions of all objects and optical 
counterparts is lower in the case of the SXDS sample, since 
the optical data are deeper in comparison with the GWS 
sample. Since X-ray data are deeper in the GWS field, and 
since it is important to deal with the complete sample when 
measuring the LF, we will use the limiting fluxes corre- 
sponding to the completeness levels of SXDS samples in 
all analyses. 

2.2 Optical data 

The SXDS optical imaging observations were carried 
out with the Suprime-Cam on the Subaru Telescope We ob- 
served five continuous sub-fields covering an area of 1.2 
deg^ with a total integration time of 133 hours (Furusawa et 
al. 2008). We downloaded and used publicly available cat- 
alogues from the SXDS web pag^!^ obtaining photometric 

http://www.naoj.org/Science/SubaruProject/SXDS/index.html 
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information in five broad-band filters B, V, Rc, i', and z, 
with limiting AB magnitudes of 28.4, 27.8, 27.7, 27.7, and 
26.6, respectively. The total number of objects in each field 
is around 900000. 

GWS observations used in this work were carried out by 
the OTELO team using the Prime Focus Imaging Platform 
(PFIP) at the William-Herschel Telescope (WhtEI cover- 
ing a total area of 0.18 deg^. Three pointings were observed, 
each of them using the B, V, R, and / broadband filters, and 
performing several exposures at each pointing of 600, 800, 
900, or 1000 sec depending on the filter Standard reduction 
procedures have been applied. Resulting limiting AB mag- 
nitudes are 24.8, 24.9, 24.6 and 23.8 in B, V, R, and / band, 
respectively. The total number of objects is around 44000. 
A detailed description of data reduction and optical source 
detection can be found in Cepa et al. (2008). 

2.3 NIR and MIR data 

NIR information for the SXDS field was taken from the 
UKIDSS DR2 catalogu|l(Hewett et al. 2006; Warren et al. 
2007; Hodgkin et al. 2009). This data has AB Hmiting mag- 
nitudes of 23.4, 23.4, and 22.9 in J, H, and K bands, re- 
spectively, and covers an area of 8.75 deg^. MIR data for the 
SXDS field was obtained with the Spitzer telescope as part 
of the XMM-LSS survey which is one of the areas surveyed 
within the Spitzer Wide-area IR Extragalactic Legacy Sur- 
vey (SWIRE) covering a total area of 9. 1 deg^. We used the 
public Infrared Array Camera (IRACp^ catalogue^ of ap- 
proximately 250700 objects, with AB limiting magnitudes 
of 21.0, 21.1, 19.5, and 19.4 in the 3.6, 4.5, 5.8, and 8.0 ^Im 
bands, respectively. 

We used publicly available NIR Ks Palomar datcfH for 
the GWS field. The catalogue contains ^ 45000 objects wi- 
th an AB limiting magnitude of 20.6 covering an area of 
0.7 deg2. Spitzer/IRAC EGS catalogufl is publicly avail- 
able and has been used in this work to gather MIR data 
for GWS X-ray emitters. The raw catalogues contain (57.4, 
48.0, 16.3, 13.6) X 10^ objects, with AB Hmiting magni- 
tudes of 23.2, 22.8, 21.3, and 21.1 in the four IRAC chan- 
nels, respectively, and a total area of 0.33 square degrees. 

2.4 Redshifts and k-corrections of optical 
counterparts 

We cross-matched the optical data with X-ray catalogu- 
es using a search radius of 3" and 2" for SXDS and GWS 
fields, respectively. Applying the de Ruiter methodology (de 
Ruiter et al 1977; Povic et al. 2009) we obtain a complete- 
ness of 99.9% and 98.3%, and a rehabihty of 76.2% and 

' ' http://www.ing.iac.es/Astronomy/telescopes/wht/ 
http://www.ukidss.org/archive/archive.html 
http://www.cfa.harvard.edu/irac/ 

http://swire.ipac.caltech.edu/swire/astronomers/data_access.html 
' ^ http ://w w w. astro. cal tech. edu/ AEGIS/ 
http://www.cfa.harvard.edu/irac/egs/ 



87.3% in the SXDS and GWS fields, respectively. Due to 
the poor reliability in the SXDS field and, in order to de- 
crease the number of spurious detections, we also applied 
the Sutherland & Saunders (1992) methodology, measuring 
the reliability for all possible matches inside the selected ra- 
dius, and finally selecting, as an optical counterpart, the ob- 
ject with the highest probability. This additional step drops 
the percent of possible fake identifications from 24% to 6%, 
in cases where multiple matches are detected with similar 
probabilities. We find 808 and 340 X-ray sources with op- 
tical counterpart in the SXDS and GWS fields, respectively. 
The final covered areas are 1.2 deg^ and 0.08 deg^ in the 
SXDS and GWS fields, respectively. These are the areas that 
will be used for solid angle estimations in Section |4] and 
cover the overlap area of all the data employed to identify 
X-ray sources and their optical counterparts, and to mea- 
sure photometric redshifts. 

We used the ZEBRA (Feldmann et al. 2006) code to 
compute photometric redshifts for optical counterparts. Due 
to the small number of previously available spectroscopic 
redshifts in the SXDS field (fifteen objects), and the lack 
of NIR data in the GWS survey (only K data available), 
photometric redshifts were also derived using the HyperZ 
code (Bolzonella et al. 2000) as an additional check. We 
implemented and used the templates from the SWIRlj*! li- 
brary (Polletta et al. 2007), including galaxies belonging 
to different morphological types, starbursts, different types 
of Seyfert galaxies, and type-1 and type-2 QSO templates. 
Comparing the results from ZEBRA and HyperZ, we ac- 
cepted only sources where the redshift results from both 
codes agreed within less than 0. 1 for z < 1 and 0.2 for z > 1 . 
Following this criterion we derived photometric redshifts 
for 308 (38%) and 203 (60%) objects in the SXDS and 
GWS fields respectively. We focus our studies on the 208 
(SXDS) and 169 (GWS) objects that have redshifts below 
1 .4 in order to derive an LF for sources with reliable red- 
shift estimates. Additionally, we compared our photomet- 
ric redshift in the SXDS field with the spectroscopic values 
recently provided by the UDS tearr0. The relation is pre- 
sented in Figure |2] (top panel). 73% and 86% of all objects 
have photometric-spectroscopic redshift differences below 
0.1 and 0.15, respectively. In the GWS field the compari- 
son between photometric and used spectroscopic values is 
presented in Figure |2] (bottom panel). In this field, 66% and 
81% of all objects show differences below 0.1 and 0.15, re- 
spectively. In the following analysis, we used the spectro- 
scopic redshifts when available, and otherwise our photo- 
metric estimates (69% and 83% for samples from the SXDS 
and GWS fields, respectively). 

The code kcorrect (Blanton et al. 2007) was used 
in both fields in order to obtain k-corrections for the op- 
tical fluxes. In contrast. X-ray fluxes were k-corrected by 
assuming a standard power law spectral energy distribution 
(SED) with r= 1.8. No intrinsic absorption correction has 

" http://www.iasf-milano.inaf.it/~polletta/templates/swire_templates.htmll 
http://www.nottingham.ac.uk/astronomy/UDS/data/dr3.html 
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SPECTROSCOPIC REDSHIFT (DEEP2) 



Fig. 2 Relation between our photometric and spectro- 
scopic redshifts in the SXDS (top) and GWS (bottom) fields. 
Solid line shows the best correlation, when zphot=zspec, 
while the dashed lines show the ±0.1 difference limits. 



been applied to the power law SED for objects considered 
to be unobscured(see Section |3]l a fixed intrinsic absorption 
of Nh = 1.0 X lO'^'^ cm^'^ was applied (e.g., Silverman et al. 
2005; Younes et al. 2012) for obscured sources. 

3 Morphological and X-ray classification 

We used SExtractor (Bertin & Arnouts 1996) and 
galSVM (Huertas-Company et al. 2008, 2009) codes to ev- 
aluate I-band morphologies for optical counterparts of X- 
ray emitters in the SXDS field. SExtractor was used for 
detection of compact sources (CLASS_STAR > 0.9), while 
we classified extended objects with galSVM. As an addi- 
tional check we made a visual classification for the largest 
and brightest objects using the IRAF/imexam(Tody 1986, 
1993) tool and isophotal contour diagrams. The final classi- 
fications include 22% compact, 30% early-type, 18% late- 
type, 23% poorly classified and 7% unidentified sources. 
Further details about morphological classifications in the 
SXDS field can be found in Povic et al. (2012). 

We obtained morphologies for X-ray emitters with op- 
tical counterparts in the GWS field using SExtractor 
(Bertin & Arnouts 1996; to select compact sources as in the 
SXDS field) and GIM2D (Simard et al. 1998) codes, and we 
performed again a visual inspection as an additional check. 
Morphology in the GWS field was determined prior to the 
galSVM code becoming publicly available. The final clas- 
sifications were carried out in the R band (the I band was af- 
fected by fringing) and included: 42% compact, 21% early- 
type galaxies, 27% late-type galaxies, and 10% unidenti- 



fied sources. The number of compact objects is significantly 
higher in this field probably due to the lower quality of opti- 
cal data compared to the SXDS field. Further details related 
with morphological classification in the GWS field can be 
found in Povic et al. (2009). 

A coarse nuclear type classification was performed in 
both fields, based on a diagnostic diagram relating the X- 
ray-to-optical (X/O) flux ratio and the hardness ratio HR 
(2-4.5keV/0.5-2keV). The X/O ratio was computed as ra- 
tio of fluxes Jo.5-4.5kcV / Fr^ where the optical flux Fji was 
derived from the SExtractor auto magnitudes in the R 
band. This flux ratio is useful to efficiently separate between 
active and non-active/Compton-thick galaxies (Alexander et 
al. 2001; Fiore et al. 2003; Civano et al. 2007), while the 
hardness ratio HR(2 — 4.5keV/0.5 — 2keV), measured as in 
equation ([T]i, is very sensitive to absorption, and thus capa- 
ble of separating X-ray type- 1 (unobscured or unabsorbed) 
from X-ray type-2 (obscured or absorbed) AGN (Mainieri 
et al. 2002; Delia Ceca et al. 2004; Perola et al. 2004; Cac- 
cianiga et al. 2004; Dwelly et al. 2005; Hasinger et al 2008). 
Full details on the methodology applied are given in Povic 
et al. (2009). Following Fiore et al. (2003), all objects hav- 
ing X/O > 0. 1 and HR(2-4.5keV/0.5-2keV) < - 0.35 have 
been classified as X-ray type-1 (unobscured), while objects 
with X/O > 0.1 and HR(2 - 4.5keV/0.5 - 2keV) >-0.35 
have been classified as X-ray type-2 (obscured) AGN. Fig- 
ure|3]shows the X-ray type classification diagram including 
the host morphological information (compact, early- and 
late-type galaxies). A large fraction of objects (59.3%) fall 
inside the region of X-ray type-1 AGN. The dashed-line 
box indicates the 'highest probability' region for finding X- 
ray type-1 AGN (Delia Ceca et al. 2004). 34.9% of our sam- 
ple (SXDS H- GWS) fan in this box while only 11% fafl in 
the lower region of the diagram (i.e. typical of stars, nor- 
mal galaxies and Compton-thick AGN). No evidence of any 
relationship between X-ray and morphological types was 
indicated since there was no clear separation of X-ray prop- 
erties based on morphological types. Nevertheless, it can be 
seen that the majority (more than 60%) of objects identified 
as compact are placed in the X-ray type-1 region. 



4 XLFs of different AGN populations 

4.1 LF measurements 

The LF is one of the most important tools for studying 
the evolution of extragalactic sources. It is defined as the 
number of objects per unit comoving volume per unit lumi- 
nosity interval, i.e. 



dVdL 



(2) 



where N is the number of objects of luminosity L found in 
the comoving volume V at redshift z. In calculating the LF, 
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Fig. 3 Relation between the X/O flux ratio and HR(2 - 
4.5keV,0.5 — 2keV), for X-ray emitters with optical coun- 
terparts in the SXDS and GWS fields. Sources belong to 
different morphological groups: compact (black crosses), 
early- (red triangles), and late-type (blue squares) galaxies. 
Solid lines separate X-ray type-1 (unobscured) and X-ray 
type-2 (obscured) regions and the area with logX/0<-l, 
where stars, non-active galaxies, and Compton thick AGN 
can be found (e.g. Fiore et al. 2003). The dashed line box 
presents the limits obtained by Delia Ceca et al. (2004), pre- 
senting the 'highest probability' region for finding X-ray 
type-1 AGN. 



we used the l/Va method (Schmidt 1968) to measure the 
space density (dN/dV) as: 



dN 1 



(3) 



with N objects having the luminosities and redshifts in the 
interval ALAz around the centre of the bin (L, z). Va{i) is 
the survey available comoving volume within which an ob- 
ject i with luminosity L{i) could be detected and remained 
in the redshift bin Az. The main advantage of the dN/dV 
estimator (compared with simply dividing the number of 
objects found by the average search volume) is that it takes 
into account the fact that more luminous objects are de- 
tectable over larger volumes than less luminous objects. 

To derive the LF ($) through the space density in the 
observed luminosity bin and surveyed solid angle (il) we 
used the following equation: 



'^(logL^z) 



47r 1 



n AlogL ^ Vail) 



N 



1 



(4) 



A bin of logarithmic luminosity has been used in this work. 
The available comoving volume Va{i) has been measured 
through the comoving distance (e.g., Hogg 2000), for the 
concordance cosmology assumed in this work. 

The statistical uncertainty of $ was estimated using the 



following expression (e.g. Marshall 1985; Boyle, Shanks & 
Peterson 1988; Page & CaiTera 2000): 



1 

AL' 



N 

E 



1 



(5) 



A Schechter (1976) function has been used to fit the LF: 



<^{L)dL = ^*ij^ 



where ^{L)dL is the number of galaxies with luminosity 
between L and L + dL per Mpc^. This approximation has 
been widely used for more than 30 year and provides a good 
fit for both local (e.g. Loveday et al. 1992; Blanton et al. 
2001) and high redshift sources (e.g. Norberg et al. 2002; 
Faber et al. 2007; Ouchi et al. 2008; Dai et al. 2009). The 
Schechter function has three free parameters which must be 
determined empirically: a, $* and L*. The quantity a de- 
fines the faint end slope of the LF, and indicates how quickly 
galaxy number counts evolve at low luminosities. It is usu- 
ally negative, implying a large number of galaxies with low 
luminosities; if a = — I the LF is said to be flat. L* is the 
so-called characteristic Schechter luminosity, and $* is the 
number of galaxies per Mpc'^ per luminosity at the charac- 
teristic luminosity. 

We computed the best fit Schechter parameters using 
a Marquardt-Levenberg algorithm (Levenberg 1944; Mar- 
quardt 1963) by means of the IDL LMFItS function, per- 
forming up to 100 iterations. Finally, the best fit parameters 
and their la uncertainty estimates, assuming Poisson statis- 
tics, were obtained. 



4.2 LF of whole AGN sample 

Our XLF analysis has been restricted to sources with 
z < 1 .4 in order to have a sample with reliable redshifts and 
an optimal number of sources in the Lx-z plane (Figure|4]l. 
Figure |4] shows this plane in the total, and soft X-rays for 
AGN in both the SXDS and GWS fields, and in the hard 
band for the AGN belonging to the SXDS field only (since 
hard range was not selected in the GWS field). In both the 
soft and total bands, the objects with lowest X-ray lumi- 
nosities belong to the GWS field, since the GWS Chan- 
dra data are deeper in comparison with the SXDS XMM- 
Newton data (see Section IZTl i. The diagrams show the X- 
ray luminosity as a function of redshift for the entire source 
sample and for different morphological and X-ray types 
(from left to right). Luminosity limits in all three bands 
for both SXDS and GWS samples are also presented. In 
all bands we can clearly see the rapid increase of luminos- 
ity at lower redshifts, and flatter luminosity distribution at 
higher redshifts. Also, we can see that the most luminous 
objects are the most distant. At 0.7 <Z2 < 1.4, we found 
more compact sources in the total and soft bands than in 
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Fig.4 From top to bottom: Relation between the X-ray luminosity in the 0.5-7.0keV, 0.5-2.0keV, and 2.0- lO.OkeV 
bands and redshift for a sample of X-ray emitters with optical counterparts in the SXDS (black stars) and GWS (red 
diamonds) fields. From left to right: This relation is shown for the whole sample of objects, and for different morphological 
(compact, early- and late-type galaxies) and X-ray types (X-ray type-1 and type-2 sources). Vertical thin, dotted line 
defines the limit between two redshift intervals that we used to study the evolution in Sections [4.21 and 14.51 Black dashed 
and red dash-dot-dashed lines present the luminosity limits of the SXDS and GWS samples, respectively. 



the hard X-rays. Early-type objects dominate late-types at 
lower redshifts in the total and soft bands, but have similar 
populations in the hard band. 

XLFs derived in the 0.5 - 7.0 keV (TXLF), 0.5 - 2.0 keV 
(SXLF), and 2.0- lO.OkeV (HXLF) bands are presented in 
Figure |5] (left panel, from top to bottom, respectively). The 
LFs were analysed in two volumes with similar numbers 
of objects in all bands. The first volume corresponds to the 
redshift range < Zi <0.7 that represents about 1/4 of the 
total volume analysed (« 0.009 Gpc^ for a covered solid 
angle), while the second volume corresponds to the redshift 
range 0.7 < Z2 <1.4 occupying the rest of the volume. The 
TXLF and SXLF were derived from both SXDS and GWS 
samples to flux limits of 9.0 x 10^^^ and 2.25 x 10~^^ erg 
cm~^ s~^ in each band, respectively. Since the X-ray GWS 
data are deeper than the SXDS data (see Section I2.lt . the 
used limiting fluxes correspond to the completeness levels 
of the SXDS samples. The HXLF has been derived from the 
SXDS sample to a flux limit of 4.0 x 10^^^ ergcm^^ s^^. 
The distributions of objects belonging to each observed lu- 
minosity bin in the total, soft, and hard bands are shown in 
Figure |5] (right panel, from top to bottom, respectively). 



In order to estimate the error in the derivation of the LF 
due to cosmic variance when both the SXDS and GWS sam- 
ples are analysed, the Cosmic Variance CalculatoJ^ (Trenti 
& Stiavelli 2008) was used. Using both the Press-Schech- 
ter (Press & Schechter 1974) and the Sheth-Tormen (Sh-eth 
& Tormen 1999) models, the error in the LF due to cosmic 
variance is smaller than 10%. 

A good fit to the LF was obtained by using a Schechter 
function. The best fit Schechter parameters are given in Ta- 
ble [T] The shape of the LF remains similar in the three X- 
ray bands (in both redshift intervals). The faint-end slopes 
(a) of the luminosity function does not show a significant 
difference with redshift. While the SXLF shows similar be- 
haviour in Zi and Z2 bins, it shows a shallower faint-end 
slope a^-l.O, suggesting slower evolution than the other 
two bands. The plots of Figure |5] taken together show that 
hard X-ray sources have higher luminosities than those se- 
lected in the 0.5-7.0 keV band, a contrast even stronger in 
comparison with the soft X-ray objects. This difference is 
also found in the luminosity density with soft X-ray sources 
showing the lowest values. 

http://casa.colorado.edu/~trenti/CosmicVariance.html 



© ?? WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



www.an-joumal.org 



Astron. Nachr. / AN (??) 



7 



' -6 
42.0 





n c "7 n i/r^\/ 
U.O-/.U Kev 

1 




0.5-2.0 keV 1 














r 2.0-10.0 keV 

: X o<z<o.7 

r X 0.7<z<1.4 





42.5 43.0 43.5 44.0 44.5 45.0 45.5 

log L[erg s"'] 



10 

6 ^ 







--i 0.5-7.0 keV J 








r 




^ — H 






0.5-2.0 keV H 

! 








L — 0<z1 <0.7 


- - - ^ 2.0-10.0 keV^ 


; - - 0.7 < 


z2 < 1.4 











Fig. 5 (Left, from top to bottom): TXLF, SXLF, and HXLF in the Zi < 0.7 (black solid line, x symbols) and 0.7 < Za < 1 .4 
(red dashed line, stars) intervals. The best-fit Schechter parameters are shown in Table [T] The arrows represent 1 a errors in 
those luminosity bins where we have data only above or below the central luminosity. In all following figures, the arrows 
have the same significance. (Right, from top to bottom): Luminosity distribution in each redshift bin in the total, soft, and 
hard bands. 



Table 1 The best-fit Schechter parameters of the XLF for a whole sample of AGN (Figure|5]l 



" Redshift range 
Number of objects 



Energy range 




Num'' 


log*!?* 


a 


logL* 


(keV) 






(Mpc~^) 




(ergs^^) 


0.5-7.0 


0-0.7 


127 


-4.49 ± 0.66 


-1.35 ±0.57 


43.69 ±0.51 




0.7-1.4 


126 


-4.76 ±0.80 


-1.85 ±0.24 


44.82 ± 0.59 


0.5-2.0 


0-0.7 


111 


-4.31 ±0.31 


-0.99 ± 0.44 


43.03 ±0.26 




0.7-1.4 


99 


-4.09 ±0.22 


-1.15 ±0.20 


43.98 ±0.18 


2.0-10.0 


0-0.7 


97 


-4.55 ±0.42 


-1.35 ±0.38 


44.25 ±0.30 




0.7-1.4 


87 


-4.32 ± 1.16 


-1.92 ±0.69 


44.98 ±0.63 



Analysing differences in evolution between Zi and Z2 
ranges we find: 

* Luminosity Evolution 

Strong XLF was found in all bands as shown in Figure|5] 
Sources typically show lower luminosities in Zi than in 
Z2 band. The Schechter fits in the three X-ray energy 
ranges lead to the same conclusion: L* shows evolution 
with redshift, having lower values in Zi interval. 

* Density Evolution 

In general, there is no evidence for any significant AGN 
density evolution in any considered band. The sources 
become brighter as redshift increases, but given the un- 
certainties, their characteristic number density <i>* asso- 
ciated to a characteristic luminosity (L*) is consistent 
between the two analysed redshift intervals. 

Strong AGN luminosity evolution was also found in pre- 
vious studies (e.g. Miyaji et al. 2001; Ueda et al. 2003; La 
Franca et al. 2005; Hasinger et al. 2005; Silverman et al. 
2008a; Ebrero et al. 2009; Yencho et al. 2009; Aird et al. 



2010 and references therein). This evolution shows that a 
larger population of high-luminosity AGN is found at higher 
redshifts, while low-luminosity AGN are dominant at lower 
redshifts. This could be a result of: (i) selection effects (e.g., 
as redshift increases, first, the detectability of low-lumino- 
sity AGN decreases, and second, the observed volume in- 
creases, thus increasing the probability of detecting a lar-ger 
number of high-luminosity sources), (ii) true evolutionary 
effects, or (iii) a mixture of both. Selection effects will be 
tested in future, using a much larger sample of active galax- 
ies, however the most luminous sources we observe at high 
redshifts we do not observe in the local universe. 

4.3 LF of compact, early- and late-type active galaxies 

We also studied the AGN XLF in relation to morphol- 
ogy. The top panel in Figure |6] shows the XLF as a function 
of morphology, while the bottom panel shows the luminos- 
ity distribution for each morphological class. The LFs of 
compact, early- and late-type objects were derived over the 
full 0.5-7.0keV range. Best-fit Schechter parameters are 
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Fig. 6 (Top): TXLF derived for different morphological 
types: compact (black crosses), early- (red stars) and late- 
type (blue diamonds) sources in z < 1 .4 redshift interval. 
(Bottom): Distribution of objects in each luminosity bin of 
compact (black line), early- (red line) and late-type (blue 
line) objects. The best-fit Schechter parameters are shown 
in Table |2] 

given in Table |2l Due to the small number of objects in 
each morphological sample the LFs were derived without 
dividing into different redshift bins. Therefore, the evolu- 
tion rates were not studied as a function of morphology in 
the framework of this paper. No significant differences were 
found within the uncertainties between LFs of compact and 
early-type objects in the total band. They have similar steep 
faint-end slopes and similar L* and $* best-fit parameters. 
The situation may reverse for late-type objects showing a 
shallower faint-end slope, suggesting slower evolution, and 
higher break luminosity. However, a larger sample of late- 
type objects is necessary to confirm the results obtained. 

These results are, in general, consistent with those of 
Georgakakis et al. (2009) who studied the contribution of 
different morphological types to the TXLF at redshifts 0.5 - 

1 .3 for the Extended Groth Strip field (EGS). Both this work 
and Georgakakis et al. (2009) found that relative to late- 
type objects, early-type hosts make a major contribution to 
the TXLF over the entire luminosity range, while at higher 
luminosities compact objects become more dominant. 

4.4 LF of X-ray type-1 and type-2 AGN 

We derive the XLF in relation to the X-ray obscuration. 
Figure Ul (left panel) shows the TXLF, SXLF, and HXLF 
(from top to bottom) of X-ray type-1 and type-2 sources. 
The numbers of sources in each luminosity bin of both types 
are shown for each energy range (right panel). Table |3]sum- 
marises the derived best-fit Schechter parameters. Due to 
the small number of type-2 objects, the LFs of type-1 and 



type-2 sources were compared at all redshifts z < 1 .4, with- 
out studying their evolution rates. Furthermore, in order to 
provide a good Schechter fit of the X-ray type-2 objects 
in the soft band, the value of the faint-end slope parameter 
was fixed to -1.5. Different values were tested, from -1.0 to 
-2.0 with a bin of 0.25, and the value of -1.5 was found to 
give the best fit. In the other two energy ranges the faint-end 
slope parameter stays free for both type-1 and type-2 ob- 
jects. 

The luminosity distributions of X-ray type-1 and type-2 
objects appear to be similar at both total and hard bands, but 
significantly different in the soft band, as expected since soft 
X-rays are especially sensitive to absorption by the central 
torus. 

In general our results suggest that type-1 sources are 
more numerous at higher luminosities, and type-2 objects 
dominate the AGN population at lower luminosities. This is 
known as the 'Steffen effect' (Steffen et al. 2003), and was 
observed in many later works (e.g., Ueda et al. 2003; Barger 
et al. 2005; Hasinger et al. 2008). Taking into account un- 
certainties, the faint-end slopes of X-ray type-1 and type-2 
sources in the 0.5 - 7.0keV and 2.0 - lO.OkeV bands appear 
to be similar at the observed redshifts. Moreover, observing 
the best Schechter-fit parameters, in general type-1 sources 
show higher characteristic luminosities (L*) then the type-2 
objects. In contrast, the values of the number density at 
characteristic luminosity do not show significant differences 
between type-1 and type-2 sources in any of the three tested 
X-ray bands. 

4.5 Evolution of X-ray type-1 AGN 

We studied the evolution of X-ray type-1 AGN. Fig- 
ure [S] shows their TXLF (top panel), as well as the number 
of sources in each luminosity bin (bottom panel). The evo- 
lution rate was studied again in Zi and Z2 intervals. Table |4] 
summarises the best-fit Schechter parameters. 

Similarly to the entire AGN sample, evidence for a clear 
luminosity evolution was found with the most luminous ob- 
jects residing at higher redshifts. Moreover, we observe a 
weak evolution in the number densities (<&*), with high- 
redshift X-ray type-1 sources having higher <i>* values. 

5 Summary 

We present the XLF for X-ray detected AGN in the 
SXDS and GWS fields. We studied the entire sample of ac- 
tive galaxies as well as different morphological and X-ray 
subsamples. In the future we will use a much larger sample 
of X-ray detected AGN from the COSMOS field to derive 
the LF in relation to the morphology of host galaxies and 
to study the evolution of different triggering mechanisms 
(work in progress). 

LFs of our entire AGN sample were derived in total, 
soft, and hard X-ray bands. We studied objects with z < 1 .4 
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Table 2 The best- fit Schechter parameters of the TXLF as a function of morphology at z < 1.4 (Figure |6]l 



" Number of objects 



Sample 



Num" 



log$* 
(Mpc^^) 



logL* 

(ergs^^) 



Compact I 79 -5.98 ±0.84 -1.72 ±0.29 44.59 ±0.64 
Early-type | 79 -6.20± 1.37 -1.98±0.31 44.46±0.97 
Late-type I 29 -6.63 ±3.38 -1.65 ±0.57 45.04 ±3.68 
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Fig. 7 (Left, from top to bottom): TXLF, SXLF, and HXLF of X-ray type-1 (blue triangles) and type-2 (red squares) 
objects at redshifts z < 1.4. The best-fit Schechter parameters are shown in Table|3] (Right, from top to bottom): Luminosity 
distribution in the total, soft, and hard bands of type-1 (blue doted line) and type-2 (red dashed line) sources. 

Table 3 The best-fit Schechter parameters of X-ray type- 1 and type-2 sources at z < 1.4 (Figure IT) 



En. range 
(keV) 


Sample 


Num" 


log<l?* 
(Mpc^^) 


a 


logL* 
(ergs^^) 


0.5-7.0 


X-ray type-1 


120 


-5.82 ±0.47 


-1.71 ±0.15 


44.69 ±0.37 




X-ray type-2 


92 


-5.29 ±0.47 


-1.56±0.32 


43.98 ±0.31 


0.5-2.0 


X-ray type-1 


115 


-5.13 ±0.31 


-1.25 ±0.3 


43.75 ±0.23 




X-ray type-2 


50 


-4.99 ±0.20 


-1.50'' 


42.68 ±0.15 


2.0-10.0 


X-ray type-1 


82 


-5.75 ±0.97 


-1.70 ±0.40 


45.01 ±0.72 




X-ray type-2 


73 


-5.05 ±0.65 


-1.46 ±0.55 


44.54 ±0.45 



" Number of objects 

^ Fixed value of a in the Schechter fit 



Table 4 The best-fit Schechter parameters of TXLF derived for X-ray type-1 sources (Figure |8]l 





Num'' 


log$* 


a 


logL* 






(Mpc-^^) 




(ergs"^) 


0-0.7 


60 


-4.74 ±0.76 


-1.37±0.67 


43.66 ±0.55 


0.7-1.4 


60 


-4.28 ±0.53 


-1.38 ±0.45 


44.50 ± 0.37 



" Redshift range 
Number of objects 
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Fig. 8 (Top): TXLF of X-ray type-1 sources derived 
in two redshift intervals: Zi<0.7 (black crosses) and 
0.7 < Z2 < 1.4 (red stars). (Bottom): Distribution of objects 
in each luminosity bin of X-ray type-1 sources in the 
Zi <0.7 (solid black line) and 0.7 <Z2 < 1.4 (dashed red 
line) redshift intervals. The best-fit Schechter parameters 
are shown in Table |4] 

(where our photometric redshifts are most reliable relative 
to available spectroscopic redshifts) inferring source evo- 
lution in two redshift intervals (0 < Zi < 0.7 and 0.7 < Z2 
<1.4) where each redshift band has a similar number of 
sources. The data were fitted with a Schechter function wh- 
ich was found to be a good approximation in all bands. 

The following results were obtained for the XLF derived 
for the entire sample: 

★ The LF shape does not change significantly between the 
two z intervals for the X-ray bands considered. 

★ Evidence for strong luminosity evolution was found in 
all X-ray bands, with higher density of more luminous 
objects at higher redshifts. 

★ Taking uncertainties into account, no significant signs 
of number density evolution was found in any of the 
tested bands. Sources become brighter as redshift in- 
creases but their characteristic number density ($*) that 
corresponds to characteristic luminosity L* appears to 
remain very much the same in both redshift intervals. 

We studied the TXLF as a function of morphology for 
sources with z < 1 .4. Three morphological types were tes- 
ted: compact, early-, and late-type with the following re- 
sults: 

★ Compared to late-type objects, early-type hosts make a 
major contribution to the TXLF over the entire luminos- 
ity range, while at higher luminosities compact objects 
become more important. 

★ No significant differences 

were found between LFs of compact and early-type ob- 
jects in the total X-ray band. 



The TXLF SXLF, and HXLF of X-ray type-1 and type- 
2 sources at z < 1 .4 were analysed, giving the following re- 
sults: 

★ The SXLF of type-1 and type-2 sources are shown to be 
significantly different due to effects of absorption of soft 
X-rays by the AGN torus. In the total and hard bands the 
XLFs remain similar. 

■k We observed the 'Steffen effect' in this work: type-1 
sources are more numerous at higher luminosities, while 
type-2 objects are more common at lower luminosities. 

★ Type-1 sources show higher characteristic luminosities 
(L*) than type-2 objects in the soft and total bands, whi- 
le at hard energies including uncertainties the properties 
stay similar. 

★ Number densities $* at characteristic luminosities do 
not show significant differences between X-ray type-1 
and type-2 objects in any of three bands. 

Finally, the evolutionary rates of X-ray type-1 AGN we- 
re studied in Zi and Z2 intervals. Deriving the LFs in each 
range, we obtained: 

★ The LFs show similar shape. 

★ A clear evolution with luminosity was obtained, with 
more luminous AGN and higher luminosity densities at 
higher redshifts. We also observed a weak evolution in 
number densities, with higher redshift sources showing 
higher $* values. 

Acknowledgements. We thank Marc Huertas-Company for a num- 
ber of valuable comments that improved the morphological clas- 
sification of analysed sources. We also thank Jack Sulentic and 
the anonymous referee for all corrections and suggestions. We ap- 
preciate all comments from the XMM-Newton team during our 
SXDS X-ray data reduction, and we also thank the SXDS team 
for making their data available to the astronomical community. 
This research has been supported by the Spanish Ministry of Econ- 
omy and Competitiveness (MINECO) under the grant AYA2011- 
29517-C03-01. MP, IM, and JM acknowledge Junta de Andaluci'a 
and MINECO through projects PO8-TIC-03531 and AYA2010- 
15169. We acknowledge support from the Faculty of the Euro- 
pean Space Astronomy Centre (ESAC). JIGS acknowledge finan- 
cial support from the MINECO under project AYA2008-0631 1- 
C02-02 and AYA201 1-29517-C03-02. JG acknowledges support 
from the MINECO through AYA2009- 10368 project. The CEFCA 
is funded by the Fondo de Inversiones de Teruel, supported by both 
the Government of Spain (50%) and the regional Government of 
Aragon (50%). This work has been partially funded by the Span- 
ish Ministerio de Ciencia e Innovaci-6n through the PNAYA, under 
grants AYA2006-14056 and th-rough the ICTS 2009-14. This re- 
search has made use of software provided by the XMM-A'ewto;i 
Science Operations Centre and Chandra X-ray Center (CXC) in 
the application packages SAS and CIAO, respectively. I RAF is 
distributed by the National Optical Astronomy Observatories, whi- 
ch are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. Virtual Observatory Tool for OPerations on 
Catalogues And Tables (TOPCAT) was used in this work. 



© ?? WILEY- VCH Verlag GmbH & Co. KGaA, Wcinheim 



www.an-joumaLor^ 



Astron. Nachr. / AN (??) 



11 



References 



Aird, J., et al., 2010, MNRAS , 401, 2531 

Aird, J., et al. 2008, MNRAS, 387, 883 

Alexander, D. M., et al. 2003, AJ, 126, 539 

Alexander, D. M., et al. 2001, AJ , 122, 2156 

Babic, A., et al. 2007, A&A, 474, 755 

Barger, A. J., et al. 2005, AJ, 129, 578 

Barger, A. J., & Hasinger G., 2005, ARA&A, 43, 827 

Bertin, E. & Amouts, S., 1996, A&AS, 117, 393 

Blanton, M. R., et al. 2007, AJ 133, 734 

Blanton, M. R., et al. 2001, AJ, 121, 2358 

Bolzonella, M., Miralles, J.-M., & Pello, R 2000, A&A, 363, 476 
Boyle, B. J., et al. 2000, MNRAS, 317, 1014 
Boyle, B. J., Shanks, T, & Peterson, B. A., 1988, MNRAS, 235, 
935 

Bmsa, M., et al. 2009, ApJ, 693, 8 
Caccianiga, A., et al. 2004, A&A, 416, 901 
Cepa, J., et al. 2008, A&A, 490, 1 

Cirasuolo, M., Magliocchetti, M., & Celotti, A. 2005, MNRAS, 
357, 1267 

Civano, R, et al. 2007, A&A, 476.1223C 

Cool, R. J., et al. 2012, ApJ, 748, 10 

Cowie, L. L., et al 2003, ApJ, 584, L57 

Croom, S. M., et al. 2004, MNRAS, 349, 1397 

Croton, D., et al. 2006, MNRAS, 365, 11 

Dai, X., et al. 2009, ApJ, 697, 506 

Davis, M., et al., 2007, ApJ , 660, 1 

Delia Ceca, R., et al. 2004, A&A, 428, 383 

De Ruiter, H. R., Arp, H. C, & Willis, A. G. 1977, A&AS, 28, 21 1 

Di Matteo T, Springel V., & Hemquist L. 2005, Nature, 433, 604 

Dwelly, T. et al. 2005, MNRAS , 360, 1426 

Georgakakis, A., et al. 2009, MNRAS , 397, 623 

Ebrero, J., et al. 2009, A&A, 493, 55 

Faber, S. M., et al. 2007, ApJ, 665, 265 

Fan, X., et al. 2004, AJ, 128, 515 

Fan, X., et al. 2001, AJ, 121, 54 

Fanidakis, N., et al. 2012, MNRAS, 419, 2797 

Feldmann, R., et al. 2006, MNRAS, 372, 565 

Fiore, R, et al. 2012, A&A, 537, A16 

Fiore, R, et al. 2003, A&A, 409, 79 

Furusawa, H., et al. 2008, ApJS, 176, 1 

Hasinger, G. et al. 2008, A&A, 490, 905 

Hasinger, G., Miyaji, T, Schmidt, M., 2005, A&A, 441, 417 

Heckman, T. M., et al. 2004, ApJ, 613, 109 

Hewett, R C., et al. 2006, MNRAS, 367, 454 

Hodgkin, S. T, et al. 2009, MNRAS, 394, 675 

Hogg, D. W. 2000, arXiv: astro-ph/9905 1 16 

Hopkins, RR, et al. 2008, ApJS, 175, 356 

Hopkins, RR, Richards, G. T., & Hernquist, L. 2007, ApJ, 654, 
731 

Hopkins, RR, et al. 2006, ApJS, 163, 1 

Huertas-Company, M. et al., 2009, A&A, 497, 743 

Huertas-Company, M. et al., 2008, A&A, 478, 971 

Hunt, M. R, et al. 2004, ApJ, 605, 625 

Jiang, L., et al. 2006, AJ, 131, 2788 

Kaiser, C. R., & Best, R N. 2007, MNRAS, 381, 1548 

Kellermann, K. I., et al. 2012, AAS, 21910906 

La Franca, R, et al. 2010, ApJ, 718, 368 

La Franca, R, et al. 2005, ApJ, 635, 864 

La Franca, R, et al. 2002, ApJ, 570, 100 

Levenberg, K. 1944, Quarterly of Applied Mathematics, 2, 164 
Loveday, J., et al. 1992, ApJ, 390, 338 
Mainieri, V. et al. 2002, A&A, 393, 425 



Manners, J. C., et al. 2003, MNRAS, 343, 293 
Marquardt, D. 1963, SIAM Journal on Applied Mathematics 11, 
431 

Marshall, H. L., et al. 1985, ApJ, 299, 109 

Matute, I., et al. 2006, A&A, 451, 443 

Merloni, A. 2004, MNRAS, 353, 1035 

Miyaji, T., Hasinger, G., & Schmidt, M. 2001, A&A, 369, 49 

Miyaji, T., et al. 2000, A&A, 353, 25 

Nandra, K., et al. 2005, MNRAS, 356, 568N 

Norberg R, et al. 2002, MNRAS, 336, 907 

Oke, J. B. & Gunn, J. E. 1983, ApJ, 266, 713 

Ouchi, M., et al. 2008, ApJS, 176, 301 

Polletta, M. et al. 2007, ApJ, 663, 81 

Page, M. J., & Carrera, R J., 2000, MNRAS, 311, 433 

Page, M. J., et al. 1997, MNRAS, 291, 324 

Perola, G. C., et al. 2004, A&A, 421, 491 

Polletta, M., et al. 2007, ApJ , 663, 81 

Povic, M. et al. 2012, A&A, 541, 118 

Povic, M. et al. 2009, ApJ , 706, 810 

Press, W. H., & Schechter, P, 1974, ApJ , 193, 437 

Richards, G. T, et al. 2006, AJ, 131, 2766 

Richards, G. T., et al. 2005, MNRAS, 360, 839 

Schechter, P 1976, ApJ, 203, 297 

Scoville, N., et al., 2007, ApJS, 172, 1 

Schmidt, M. et al., 1968, ApJ, 151, 393 

Sekiguchi, K. et al., in preparation 

Silverman, J. D., et al. 2008a, ApJ , 675, 1025 

Silverman, J. D., et al. 2008b, ApJ, 679, 118 

Silverman, J. D., et al. 2005, ApJ, 618, 123 

Simard, L. 1998, ASPC, 145, 108 

Sheth, R. K., & Tormen, G., 1999, MNRAS , 308, 119 

Steffen, A. T, et al. 2003, ApJ, 596, L23 

Sutherland, W., & Saunders, W. 1992, MNRAS , 259, 413 

Tody, D., 1993, Astronomical Data Analysis Software and Sys- 
tems II, A.S.P. Conference Ser., Vol 52, eds. R.J. Hanisch, 
R.J.V. Brissenden, & J. Barnes, 173 

Tody, D., 1986, Proc. SPIE Instrumentation in Astronomy VI, ed. 
D.L. Crawford, 627, 733 

Trenti, M., Stiavelli, M., 2008, ApJ , 676, 767 

Ueda, Y, et al., 2003, ApJ , 598, 886 

Warren, S. J., et al. 2007, MNRAS , 375, 213 

Wu, Y., et al. 201 1, ApJ, 734, 40 

Yencho, B., et al. 2009, ApJ , 698, 380 

Younes, G., et al. 2012, A&A, 539, 104 



www.an-joumal.org 



© ?? WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim 



